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ABSTRACT 



A method for resolving the mesoscale distribution of precipitation over mountainous 
regions is developed using synoptic-scale aneilyses and topography with a 40 km 
resolution. The method consists of calculating the saturated moisture flux at a level 
assruned to be tire cloud base and correlating this quemtity with observed precipitation. 
General correlation coefficients ranged from -0.2 under conditions poorly represented by 
the synoptic surface analysis to -0.55 when these analyses are more representative of actucil 
forcing. When applied, regionally, correlation improved dramatically, suggesting potential 
incorporation into precipitation forecasts over mountainous areas. Over wider regions of 
less defined, topography, correlation is weaker but is still an improvement over the current 
synoptic-scale precipitation forecast products. 
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L INTRODUCTION 



TTie enhancement of synoptic-scale precipitation by topography is well understood, 
and has been routinely verified. Fine-mesh mesoscale models have proven particularly 
accurate in reflecting orographic influences on large-scale precipitation events (Colton 
1976, Collier 1975). Simpler estimates of this enhancement factor have been shown to 
be almost exclusively functions of regional topography combined with low-level wind 
fields, synoptic-scale vertical motion, and a knowledge of the low-level relative humidity 
(Browning et al. 1975, HUl et al. 1981). 

Thus far, only detailed mesoscale models have been able to reproduce the small- 
scale orographic enhancement process. With synoptic-scale precipitation products 
showing little or no regional structure, and mesoscale forcing fields generally unavailable, 
forecasters rely mostly on local knowledge and vague statistical multipliers to predict 
areas of heavy rainfall. Consequently, the ability to quantify local orographic effects and 
theii associated distribution of precipitation based solely on synoptically derived wind and 
moisture fields would be of particular use to the forecaster. 

In this study, the above mentioned quantities of subsynoptic-scale topography and 
synoptic -scale wind and moisture data are combined in an effort to reproduce areas of 
enhanced precipitation diurmg a series of eastern Pacific extratropical cyclones which 
passed over California during the period 1200 TJTC 27 Pebruary- 1091 to 1200 UTC 04 
March 1991. Topoji'ar''^'. coti’nine'l fi- i.i i-: by the 
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National Meteorological Center’s (NMC) final analyses, produces a distribution of the 
total vertical moisture flux. This vertical moisture flux estimate is correlated with hourly 
precipitation data collected by the National Oceanic and Atmospheric Administration 
(NOAA), covering 244 reporting stations throughout California for the same period, to 
produce an equivalent rainfall amount. 

The overall intent of this study is to aid forecasters in identifying areas most likely 
to experience orographic enhancement of precipitation for a given synoptic event and 
pjovide a quantitative precipitation estimate. This procedure may be of use in estimating 
seasonal snowpack levels in otherwise inaccessible regions, as well eis providing near term 
flood potential for specific drainage basins, while only requiring synoptic-scale 
infonnation. More specific objectives of this study are to: 

* Assess the degree of accuracy to which precipitation distribution can be represented 
using mesoscale topography and large-scale wind and moisture fields. 

- Develop a simple quantity, calculated locally from model output products, to 
represent the enhancement of precipitation due to topography. 

* Present specific synoptic-scale conditions and/or regional topographic characteristics 
which either promote or discourage the use of this parameter. 

* Determine the level of improvement, if any, the employment of this parameter 
shows over current precipitation forecast products. 

Following a discussion of background work done in the field of orographic 
precipitation enhancement, an outline of the dynamical basis for the vertical moisture flux 

as a precipitation pa.ra.n’eter pte.'^en'c'!. a - -v- n < -• t (!!►' I'iquited 

input fields. Next, thi.5 para.uictet i," -vi 'il' <'n'<ot-\(fd (''‘•crnaiion to obtain a 
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statistical fit assigning specific precipitation accumulations to the calculated values 
followed by an examination of the factors which affect correlation. Results obtained by 
applying the precipitation parameter both generally and regionally are examined, including 
a comparison to current precipitation forecast aids (specifically, the Nested Grid Model 
(NGM) 12 hour precipitation forecast). Conclusions focus on the usefulness of this 
technique as a forecast aid. 
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n. BACKGROUND 



It is certainly no secret that mountainous areas often receive greater amounts of 
precipitation than do upwind or downwind low-lying regions. However, it was not until 
Bergeron (1965) introduced his "seeder cloud - feeder cloud" concept that a specific 
mechanism for topographically enhanced precipitation became accepted. As described by 
Raddatz and Khandekar (1979) this process consists of low-level (feeder) clouds forming 
in upslope areas which "wash out" when seeded by raindrops falling from synoptic -scale 
(seeder) clouds above. The idea is that the feeder cloud is constantly replenished and will 
continue to enhance precipitation in proportion to its depth, provided seeding from above 
is uninterrupted. Storebo (1976) and Bader and Roach (1977) were able to confirm 
Bergeron’s theories by results which related a 10% increase in precipitation to a 1000 m 
feeder cloud depth. Interestingly, Browning et al. (1974) suggested that in some instances 
(namely, cases of mid-level potential instability), orography alone could be responsible 
for generating precipitation, even if synoptic-scale seeding was not available. According 
to Browning et al. (1975), condensation over mountainous regions is roughly proportional 
to the low-level wind velocity component normal to the topographic barrier combined 
with synoptic-scale vertical motion and moisture quantities. This relationship was 
generally verified by the mesoscale numerical model developed by Collier (1975). 

As mentioned in the introduction, Colton (1^76) achieved retnarkably accurate 
results using a fine-mes'i (‘•'.5 kn) s:iid) :• i ;T>Mi.irtT t-- ' • i'- • 1 1 ^ 75 ). applied 
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to the Sierra Nevada of California. A key assumption made by Colton was to define the 
orographic precipitation as a quantity dependent on the presence of synoptic-scale 
precipitation. As such, his scheme is only able to enhance synoptic-scale precipitation, 
and cauinot force precipitation based solely on orography. 

A later study by Hill et al. (1981) further confirmed the seeder-feeder enhancement 
theory. They additionally showed that if the total precipitation rate exceeded a given 
threshold value (0.5 mm hr '), the topographic enhancement factor was not dependent on 
the background precipitation rate, .hr short, as long as the background rainfall rate was 
sufficiently high, enhancement was simply proportional to the low-level moist flow 
interacting with topography, and was not a function of the amount of seeding from above. 
This carries the important implication that orographically induced precipitation can simply 
be added to synoptic-scale rainfall to obtain a total, rather than the orographic term being 
a function of the synoptic -sc ale term. This relationship provides the basis for deriving 
the precipitation parameter described in the following chapter. 

More recent studies have attempted to modify feeder-seeder theory to account for 
additional effects including stratification and wind drift. While showing the effect of 
stratification to be negligible as it relates to total washout, Camithers and Choularton 
(1983) did show wind drift effects, which were not accounted for in models such as Bader 
and Roach (1977), to be detrimental to both local maxima and total enhancement by 
"advecting" orographic precipitation to downwind locations. Additionally, this study 

backed up the idea presepied bv <■/ 'rl, i !'v.> > ■!’ >• '-niinti- .. nujepeu'lent of 

seeding rate (also contradictoi'v to tho P-aHet t'tvi ( 1077 ' moii'ci' 



5 



Most recently, Stow. Bradley and Gray (1991) conducted a study with results 
largely consistent with seeder-feeder theories, and emphasized the significant differences 
between the previous enhancement studies. They regarded the work so far conducted as 
inconsistent in terms of dealing with the effects of the low-level wind speed, seeding rate, 
and barrier height. Moreover, they felt that although the basic physics of topographic 
enhancement are well understood, refinement is required to handle the above 
inconsistencies. The present study is not concerned with any of these "refinements", but 
rather with using the existing theory to derive a simplified method of representing 
mesoscal e-resolution topographically enhanced precipitation from exclusively synoptic- 
scale data. 

It should be noted that topography plays a larger role in weather than simply 
enhancement of precipitation. In fact it has been shown that topography can influence 
the low- to mid-level flow in such a manner as to completely reposition large-scale 
pressure centers (Sawyer 1959), as well as greatly altering the small-scale wind patterns 
(Danard 1977). An attempt to combine small-scale deformation of the wind field due to 
topography with established enhancement theory was undertaken by Anderson and Nilsson 
(1990) and had only moderate success, mostly owing to large errors in obtaining 
precipitation measurements m drifting snow. Any deformation of the wind field by 
topography will be ignored in this thesis, in terms of calculating our precipitation 
parameter. As such, a main objective of this study is to determine the feasibility of 
representing mesosca.ie fli-.ii n. iv,- •• ihe actual 

mesoscale wind fields 
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m PROCEDURE 



A. THEORY 

Preceding studies have estimated orographic enhancement of precipitation by 
determining the vertical extent of a feeder cloud which forms in upslope regions based 
on low-level winds and wet bulb potential temperature. This cloud depth can be directly 
related to orographicaJly induced precipitation provided seeding from above (synoptic- 
scale rain) is available. Assuming the "washout" efficiency is properly handled, numerical 
models can easily determine the enliancement factor based on calculated feeder cloud 
depth (Raddatz and Khandekar 1979, Bader and Roach 1977). 

Since the purpose of this study is to provide an easily-calculable means of 
estimating topographical enhancement of precipitation, a slightly different approach than 
a numerical model simulation is taken. This estimate of topographic precipitation 
enhancement is based on a "ground-up" quantity reflecting the total saturated moisture 
flux at a given height (assumed to be the cloud base), (wq)„. Conceptual restraints are 
imposed such that any upward saturated moisture flux at this level is condensed, and not 
allowed to either be redistributed into higher levels of the atmosphere, or result in 
supersaturation. Essentially, any such flux is completely "washed-out" as precipitation. 
Since this flux is representative of precipitation only under saturated conditions, a 

threshold 700 m.b relative' buvnidhv vai"‘- iii\ jj; to 

represent saturation at the svnoptic scalf"* ■'■'■'ll n*' ,>i'i'lb'i ib-' '"old. 
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B. DERrVATION OF PRECIFXTATION PARAMETER 



Bader and Roach (1977) used the following form of the continuity equation to 
represent the liquid water content: 



Dq^ 

Dt 



=C-A, 



( 1 ) 



where DqJDt represents rate of change of liquid water content following the air parcel 

trajectory. C is used to represent the condensation rate while A is the rate of washout 
from above. It follows that Dq JDt is simply the total precipitation minus evaporation 

under saturated conditions, such that 



Dq. 

Dt 



—{P-E), 



( 2 ) 



where P and E stand for precipitation and evaporation, respectively. 
(1979, pp 40-44), . ^ 



dq.. ^ ^q.. 

dt dt 



dq,. ^q,. 

+v 

dx dy 



dj_ 



J 



Following Holton 



represents the total time derivative of q^, in terms of the local tendency and advective 

components. Rearrangement in vector notation and use of the isobaric form of the 
continuity equation with (2), gives 
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^^ViqU)=-{P-E). 

dt 



( 3 ) 



Changes in the saturation specific humidity {qj result from temperature and pressure 
changes, which are largest due to vertical advection. Thus, the local tendency dqjdt is 

smaller than the flux divergence "V-qU and will be ignored in this calculation. The P-E 
term is the net condensation, which is equivalent to the rainfall when integrated over the 
depth of the atmosphere. The net condensation is written simply P in the following 
equations. We are left with 



-P^^iqU), 



( 4 ) 



as a first-order approximation to the net precipitation. Separating (4) into horizontal and 
vertical components yields 



-F = 



dx dy 



. djqo)) 
dp 



This can be integrated from the surface to the top of the atmosphere H, to give 
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( 5 ) 




Ph 

where (^co)^ represents total moisture flux at H, which is zero. The term - 

p„ 



represents the rainfall. R, and the right-hand side of equation (5) gives the contribution 

Ph 

from both topographic uplift, (<?co)^ , and synoptic-scale forcing, + ^\qU)]dp . The 

Po 




synoptic -scale forcing is simply the net horizontal flux divergence of saturated air into the 
column. Because is largest at low levels, the net moisture flux is dominated by the 

lower part of the atmosphere. Consequently, this term is approximated by the integrated 
moisture flux convergence from the surface to 700 mb. The approximate rainfall is then 
given by 

700 

To 

with R representing net precipitation in the colnmn noth terms on the right- 

hand side of the equation are easily obtained .^taodatci vv itvi :nvi nioisture data 
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combined with surface topography. Tlie 40 km resolution topography used in this study 



'^h 9 /} 

is shown in Figure 1. Specifically, the term (^co) is given as q [m__+v — ], 

dx dv 



with u, V and h representing horizontal wind velocity components and surface slope, 
respectively, while q^ refers to the surface moisture value. 

The surface forcing (^vr) , and ultimately {q(0) , was calculated using the 1000 mb 

/’p Po 

wind field and topography in meters, with the 1000 mb winds treated as crudely 
representative of synoptic-scale surface winds for this purpose. The resultant surface- 
forced moisture flu:x is in units of gfkg-cmfs. An approximate conversion to pressure 
coordinates is made by applying the relation qw=-pg(x)q . 

700 

The synoptic scale term ( + is calculated using a centered finite 

r„ 



difference scheme to compute the moisture flux divergence at levels specified by the wind 
and moisture field inputs. Summations over the desired vertical distance yield the net 
moisture flux divergence in the surface to 700 mb layer, which is the vertical moisture 
flux q(o through the 700 mb level. This is the quantity wOf, referred to in the Results 

chapter. Figures 2 and 3 show sample contoured surface vertical velocity and qa> flux 
fields obtained frorn these methods f---' "le — -'i- ' in Ficnie 2 has 

been subjected to a five passes tlv'^ngh a fine'"-p''"'' ;'”oo'hi.nc i”'-'.- ''; to ensure this 
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quantity is of the same order of magnitude as the synoptically derived forcing in Figure 
3. Finally, the total moisture flux at 700 mb, (hereafter referred to as wqf), is 

obtained by simply adding the topographic and synoptic-scale components, utilizing the 
rough equivalency of pbars ' and cms ' as described above (Figure 4). Thus both wO/and 

wqfaie in units of gkg "1 fibars ’I . This unit definition applies to all figures showing wOf 
or wqf fields unless otherwise noted. 

Note that this {qto)^^ quantity is not in units of rainfall. The actual associated 
rainfall is obtained by first employing a 75% 700 mb relative humidity threshold to the{<x>q)^^^ 

field to ensure saturation. This field is then correlated with observed 12 hour 
accumulated precipitation centered on the analysis time. A linear regression is cdculated, 

yielding a conversion factor based on the slope and y-intercept, allowing 

quantities to represent 12 hour precipitation totals. For the overall correlation, no 
objective limitations are imposed other than the relative humidity threshold. Regional 
constraints are imposed, however, when examining the effectiveness of this technique in 
more localized mountainous regions. 

Perhaps the most obvious limitation of using this quantity to estimate precipitation 
is its exclusion of sub-grid scale (i.e., convective) precipitation processes. While no effort 
is made to quantify the effect of topography on such precipitation, these conditions are 
identified using satellite imagerv'. a.uH o.''e noifH .n<; in .....Uv.:;.; pnq results 

section. 
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C. DATA 



The saturated moisture flux quantity used to represent precipitation is calculated 
from NMC final Global Data Assimilation System (GDAS) analysis fields, which were 
then statistically fitted to hourly precipiiaiion ru’''!-'- moaa. The 
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NMC GDAS analyses include twelve pressure levels from 1000 to 50 mb (as described 
by Kanamitsu 1989). For this study, the 2.5" NMC GDAS grid is interpolated to the 40 
km grid on which the regional topography was derived. The wind and moisture fields 

at the 1000. 850 and 70 n mi' l^'v^'i': ' '<• !i>- ' ’•"lyied n<oi.'^ture 

flux. H-qf. 
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Figure 3. SynopticaJJy forced 700 mb moisture flux (wO/). 1200 UTC 01 March 
1991. Negative values correspond with upward flux. 
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Figure 4. Total 700 mb moisture flux {wqf), 1200 UTC 01 March 1991. 
Negative values correspond with upward flux. 
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Specific rainfall totals obtained through hourly precipitation data publications 
(NOAA 1991) for the period of study were used to correlate precipitation with the wqf 
values. NOAA’s precipitation totals consist of hourly rain gauge observations made at 
244 locations within California. For our purposes, the observations have been grouped 
in six-hour increments exactly preceding and following the NMC analysis time. This 
gives 12 hour precipitation totals centered on the NMC analysis time which are then 
correlated with the corresponding values. 
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XV. SYNOPTIC ANALYSIS 



A. OVERVIEW 

The series of synoptic-scale events that occurred during the period of this study 
were produced by the development of a 500 mb trough off the Northern California Coast. 
The system center migrated northeastward, over the period, and supported a series of 
short-wave troughs and associated surface cyclones (Figures 5 and 6). 

An important factor in the synoptic evolution of this period is the movement of the 
primary jet stream. Approximately midway into the period (1200 UTC 02 March 1991), 
the jet repositioned itself over Central California, from the upper Baja Peninsula where 
it was located early in the period. This shift is evident in the 300 mb isotach analyses 
representative of the jet position at 1200 UTC 28 February 1991 and 0000 UTC 03 March 
1991 (Figure 7). The jet axis shifts from a westerly to a more southwesterly orientation 
coincident with its latitudinal shift. As a result, the weather affecting California was 
dominated by highly baroclinic, fast moving cyclones in zonal flow with the surface 
center well poleward of the maijn jet early in the period, and weaker, less baroclinic 
systems late in the period, with only a brief transitional time in between. 

As a consequence of the poleward migration of the jet, incoming weather systems 
affected more Northern sections of California late in the period, as evidenced by the 500 
mb height rises obse'^ed ' (') ilie 
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raix) is concentrated m Southem Califoniia early in the period, with virtually no rain 
observed in Southem California by the end of the period. 
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Figure 5. 500 mb heights (meters, left column) and surface pressure, including 
frontal positions (mb, right column) for (from top) 1200 UTC 27 Febmary 1991, 
1200 UTC 28 February 1991 and 1200 UTC 01 March 1991. 
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Figure 6. 500 mb heights (meters, left column) and surface pressure, including 
frontal positions (mb, right column) for (from top) 1200 UTC 02 March 1991, 1200 
me 03 March 1991 and 1200 l.rrc 04 March 1991 
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B. CVCLONE CHARACTERISTICS 



I. 0000 UTC 28 February 1991 - 1200 UTC 02 March 1991 

As mentioned, cyclones affecting the weather in California during this period 
developed on the cold, cyclonic side of the main jet, as seen in the GOES infrared 
imagery at 1200 UTC 28 February 1991 (Figure 8) over the Eastern Pacific and the 
conesponding mean sea level pressure and 1000-500 mb thickness analyses (Figure 9). 
This particular image and analysis are indicative of the general nature of the cyclones 
during this period. The strong surface thickness gradient along 30° N as well as the low- 
level instability upstream of the main frontal cloud bands on the GOES imagery indicate 
the strength of the baroclinic systems passing through California. These systems had 
significant pre-frontal cloud and precipitation bands associated with them. In general, this 
tyi?e of pattern is similar to that in which Collier (1975) applied his mesoscale 
enliancement model over North Wales, thus providing a good basis for comparison. 

Cyclones of this nature are favorable for observing orographically enhanced 
precipitation, largely due to the forcing being well represented by the large-scale analyses 
and the precipitation primarily associated with the occluded front. However, post-frontal 
convection does pose a problem. While CoUier (1975) applied his technique to rain 
associated with both cold and waim frontal passages, Bader and Roach (1977) tried to 
focus only on the mostly stable precipitation produced in the warm sector of cyclones. 
Browning et al. ( 1974 ) earlier concluded that orograrhicaUy induced rrecipitation was 
most evident in the v a^n a." phu-f i > i:-,. .. ,,,|p duiing tutd 
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after cold frontal passage. Still, tbek results indicated orographic effects could be 
quantified even during and after the passage of the cold front. Results obtained by Collier 
(1975) showed less degradation of quantifiable orographic enhancement during cold 
frontal passage, as compared to warm sector precipitation, than did the results presented 
by Browning et al. (1974). Collier (1975) placed more emphasis on distinguishing 
between fast moving baroclinic systems in which mesoscale forcing is rniiiimized, and 
slower, less baroclinic features as fat: as the observability of orographic precipitation 
enhancement. Since the cyclones in this period of the present study can be characterized 
as occluded surface features with mostly stable precipitation prior to frontal passage, and 
some post-frontal convection, orographic effects should be well represented. This study 
does not examine in detail the effect of mesoscale post-frontal convection on the overall 
correlation between wgf and observed precipitation. 
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Figure 7. 300 mb heights (solid, meters) and isotachs (dashed, knots) for 1200 UTC 
28 Febniaiy 1991 (top) and 1200 UTC 02 March 1991 (bottom). 
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2. 1200 UTC 02 March 199.1 - 1200 UTC 04 March .1991 



The general nature of the three weather systems affecting California during this 
period is represented in Figures .10 and 11. The thickness gradient has shifted northward 
off the California coast and the amplitude of the pre-frontal thickness ridge is less than 
in the earlier cyclones. In addition there is a lack of post-frontal instability depicted on 

die infirared image. An exception to this pattern is the area of cloudiness at approximately 35° 

N, 150° W. As with the systems preceding it, this particular system becomes cut-off 
from the cold air within which it developed, and evolves mostly in response to the strong 
southwesterly flow aloft. This system (which reaches California at approximately 1200 
UTC 04 March 1991). like the two seen ahead of it, has little representation of a 
thickness ridge ahead of the front as it reaches land. Without the pre-frontal baroclinic 
forcing, the rainfall associated with this type of system is less dominated by stable pre- 
frontal precipitation than cyclones in the earlier period. Due to the uninterrupted supply 
of warm, moist low- and mid-level flow from the southwest during this period, incoming 
weather systems produced less organized and more convective -type precipitation. This 
type of pattern was specifically avoided in terms of quantifying orographically induced 
precipitation by Collier (1975), since the forcing responsible for the distribution of 
precipitation was difficult to represent, even on the mesoscale. Consequently, this period 
of the present study may not be adequately represented by the synoptic-scale and 
orographic effects on precipitatiou. 
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Figure d. GOES infrared imagery. 1200 UTC 28 February 1991. 
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3- Summary 



The above synoptic discussion suggests two cases that produced the observed 
rainfall distributions: 



* Case (i) - strongly baroclinic cyclones with a well defined pre-frontal thickness 
ridge. 

- Case (ii) - subtropically influenced convective complexes with little or no pre- 
frontal thickness ridge. 

The first case tends to dominate the winter rainy season over California, 
although, the second case is not uncommon. Consequently, results from the two periods 
in this study cover contrasting synoptic conditions which produce rain in California. 
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V. RESULTS 



A. OVERALL 
X- Statistics 

Table 1 shows the correlation, slope, and y-intercept from a linear regression 
analysis of wqf znd observed 12 hour precipitation for four representative analysis times. 
The correlation values run from a virtually non-correlative 

-0.19 to a marginally respectable -0.55, as shown by the representative scatter plots in 
Figures 12 and 13. The correlation values are grouped at -0.2 and -0.5 for the later and 
earlier synoptic periods, respectively. This suggests that the synoptic and mesoscale 
factors discussed in the previous chapter predetermine the success of this technique to 
be either moderately correlative (-0.5) or essentially non-correlative (-0.2), with little in 
between. Averaging the slope and y-intercept of all "correlative" times (Table 1 shows 
just representative times) yields -0.010 and +35 respectively. Therefore, 

[-(m'^x0.01]+35=F 

equates calculated wqf values with expected 12 hour accumulated precipitation (P), in 
inches. Note that all values of P less than zero represent a precipitation value of zero. 

As will be shown in a later section, correlation is greatly improved in 
mountainous subregions as comp-aved to the ease, tech.nique is 

intended to be "univetsal’’ in its opi'lica'inn ei'' •• :calc parameters 
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and topography, (6) will be applied in aU sub-regions, rather than deriving specific 
statistics for individual sub-regions. 

TABLE 1. LINEAJR REGRESSION STATISTICS {WQF VS. 12 HOUR 
PRECIRTTATION) 





CORRELATION 

COEFFICIENT 


SLOPE 


Y- 

INTERCEPT 

(INCHES) 


0000 UTC 28 
1 FEB 1991 


-.55 


-.013 


+0.50 


1200 UTC 01 
MARCH 1991 


-.54 


-.010 


+0.35 


0000 UTC 03 
, MARCH 1991 


-.25 


-.003 


+0.40 


1200 UTC 04 
MARCH 1991 


-.19 


-.003 


+1.00 



One interesting statistical result is the fairly consistent y-intercept value for the 
linear regression of about +35 (average of six "correlative" times). This indicates that the 
precipitation threshold (the x-intercept) generally corresponds with sUghtly downward 
saturated moisture flux values. Though the smoothing process tends to distribute upward 
flux values into lee slope regions, it is still likely out- weighed by the "spillover" effect 
observed by Colton (1975). In this, spillover process, upward vertical velocities are 
advected over the summit region, effectively "mixing" upward vertical velocity air into 
lee regions. This process is not accounted for in the wqf calculation, which sees all lee 
regions as having downward ve1r>chies in pS'r.'tx >i. "pn ard 'noiinti. 
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Figure 9. Surface pressure (mb, solid) and 1000-500 mb thickness (meters, 
dashed). 1200 UTC 28 Februaiy 1991- 

AdditionaLly, horizontal wind drift, as described by Collier (1975), may be producing 
artificially high precipitation 

observations across summit regions, especially under conditions of strong low-level flow 
with a substantial perpendicular-to-barrier component. 

The other consistent statistical feanire is the slope value of either about -0.01 
for the correlated times, or about -0.003 for times showing less correlation. This suggests 
a "flatter" precipitation distribution during non-correlative times, with non-orographic 
pTocesses rnore ^7^ frM< 
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2. Distribution 

Although the lineaj regression statistics computed above have a relatively low 
correlation, comparison of horizontal distributions of the w^-inferred precipitation and 
observed precipitation reveals qualitatively correct distribution. In fact, this comparison 
reveals a significant improvement over NGM model forecast precipitation. 

a. 0000 UTC 2S February 1991 - 1200 UTC 02 March 1991 

Figures 14 and 15 depict the statewide precipitation distributions produced 
by M'0/(the synoptic-scale forcing described in Chapter HI), NMC’s Nested Grid Model 
(NGM) forecast precipitation for 1200 UTC 28 February 1991, actual observations and 
wqf. Of note here is the similarity between the NGM forecast precipitation and the wOf 
field, particularly the maximum over Southern California. This implies that the synoptic 
forcing we are using in our calculations of wqf bears rough equivalency to the NGM 
precipitation. This time is especially informative since the 12 hour 700 mb vertical 
velocity forecast verifies paulicularly well compared to the kinematically calculated 
analyzed 700 mb vertical velocity (not shown). It should be mentioned that wOf showed 
virtually no correlation with observed precipitation throughout the study. 

In general, we see a significantly greater amount of detail within the 
precipitation field produced by wqf, compared to either uO/ or the NGM forecast. Also 
note the "false" maxima produced in the observed field. This is mainly a result of the 
multiquadrant contouring technique extrapolating basis functions into regions of sparse 
data (i.e., Mojave .Oesevt, ove' the ^ h- :■ ••• .i-; :> ;id-- ii>. o( the 
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Figure 10. GOES infrared imagery 1200 UTC 02 March 1991. 

distribution reflects the observations reasonably well, particularly the maxima observed 
over the mountains of Southern California and the Sierra Nevada. Note that correlations 
are performed between the gridded wqf fields and station locations (rather than gridded 
precipitation). It should be noted that overall correlation is generally constant from 0000 
UTC 28 February 1991 to 0000 UTC 01 March 1991, with the exception of time 0000 
UTC 01 March 1991, which has a correlation coefficient of -0.21. This drop in 
correlation is most likely in response to the low-level convection seen in Figure 8 located 
between the main bn'Kl". ’T'lvs pnM''”' •• > : - ] "tndies. that 
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post-frontal precipitation distribution is determined less by orography than convective 
processes. 

b. 1200 UTC 02 March 1991 - 1200 04 March 1991 

Figures 16 and 17 depict the same parameters as Figures 14 and 15 
except for 1200 UTC 04 March 1991, vvhich is representative of the second period. Here 
H'^is generally a less accurate representation of observed precipitation, particularly along 
the Central California Coast. Still, orographic maxima produced by wqf over the 
mountains of Northern California and the Northern Sierra Nevada are reasonably 
representative of observed maxima. Though the corresponding NGM precipitation 
forecast shows good areal coverage, it is devoid of any of the above structure. In this 
case the NGM vertical velocity forecast, on which precipitation is largely based, shows 
a strong area of upward motion located over Northern California, whereas the verifying 
kinematically calculated analysis shows this area of upward motion to remain off-shore 
(not shown). As such, the NGM precipitation field owes its accuracy more to chance than 
to improved dynamical representation. It seems that neither the NGM forecast or the 
GDAS analysis fields are able to pick up on the forcing responsible for the locally heavy 
rainfall in tliis period. Referring back to Figure 6, we see em analyzed "front" off the 
Notthem California coast at this time, with most of the state within the warm sector. A 
random check of station observations supports the strong southerly flow depicted by the 
GDAS analyses. Convective lifting would, in this case, be the best explanation for the 
lack of upward vertic.?! veto'.uv vn \ A'; noif'l 
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Figure 11. Surface pressure (mb, solid) and 1000-500 mb thickness (meters, 
dashed). 1200 UTC 02 March 1991. 



earlier, this period is dominated by only marginally baroclinic features in the moist 
southwesterly flow aloft. As such, convection is likely to be equally important in both 
the pre- and post-frontal periods. Though convective processes are parameterized in the 
NGM precipitation forecasts, they are not accounted for at all in the calculation of wqf. 
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3 - Correlation. Lirnitiog Factors 



o. Synoptic Conditions 

As described by Collier (1975), synoptic conditions favorable for 
observing orographic enhancement of precipitation are characterized by fast moving 
baroclinic systems. The thinking is that slower moving, or non-baroclinic systems allow 
more non-orographic and/or mesoscale mechanisms to influence rainfall distribution. This 
is not to say that topography does not influence mesoscale structure. Rather, since only 
the large-scale structure is considered, topography is only allowed to influence the 
synoptic-scale fields. As mentioned previously, the period prior to 1200 UTC 02 March 
1991 was characterized by baroclinic features well represented by the synoptic fields, 
whereas non-orographic and/or a mesoscale processes dominated after this time. The 
statistical correlations support this characterization of the dominant precipitation type 
during the periods. However, as seen in the horizontal distribution of wqf during the 
second period, the qualitative distribution of precipitation is realistic in the second period. 
This suggests that the amount of precipitation is poorly represented in more convective- 
type periods, but the general distribution is still strongly influenced by the topography. 

b. Smoothing 

Aside from subgrid-scale effects, smoothing of the topographic w used 
to calculate ( 0 x 7 ) seems to be most detrimental to correlation values. As previously 

P 0 

discussed, iTiesoscalc’ •' 'ptic-^cp^e inoo 

nib wuids over ^0 .tu'’ rcscKnipn tppp£;ipp^iv b* * -pps'^ (o be 
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Figur€ 12. Observed precipitation (vertical axis, inches) vs. wqf (horizontal axis) 
and calculated linear regression for 0000 UTC 28 February 1991 (top) and 1200 
UTC 01 March 1991 (bottotn'). 
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compatible with the synoptic-scale forcing. As a result, stations upwind of topographic 
barriers, such as Fresno (station FRES, Figure 18) are generally assigned values far 
in excess of the actual calculated forcing. Note (Figure 18) that the smoothed wqf value 
for FRES (-15) is nearly twice that of the unsmoothed value (-8). 
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c. Location 



Certain locations, such as coastal areas, are more likely to have higher 
actual relative humidities than inland areas, which are analyzed as equal on the synoptic- 
scale. As a result, more of the vertical moisture flux wqf is likely to be converted to 
precipitation in coastal areas than inland areas, where a larger amount of moisture flux 
{wqf) may be required to maintain or reach saturation. This relationship is evident in 
comparing the linear regression corresponding to the San Bernardino - San Gabriel 
subregion (subregions are defined in the following sections) to that of the Southern Sierra 
Nevada (Figure 19). Note that a value of -100 represents approximately a 2.75" 12 
hour precipitation accumulation according to the linear regression derived for the San 
Bemardino/San Gabriel subregion, but less than 1.0" based on the Southern Sierra Nevada 
subregion’s statistics. 

d. Other Considerations 

Obviously, numerous other considerations affect the correlation of 
precipitation and 700 mb saturated moisture flux. Redistribution of moisture above 700 
mb, quantification of evaporation, and wind channeling induced by the topography are a 
few examples. Though important, such processes cannot be easily accounted for on the 
synoptic scale. In fact, these factors will require more specific consideration if wqf 
cannot adequately represent the forcing patterns which lead to precipitation distributions 
in mountainous regions. 
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Figure 1,3. Observed precipitation (vertical axis, inches) vs. wqf (horizontal axis) 
ajid calculated linear regression for 0000 UTC 03 March 1991 (top) and 1200 UTC 
04 March 1991 (bottom). 
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Figwrt 14. NGM 12 hour precipitation forecast, valid 1200 UTC 28 February 
1991 (mm, top) and uO/ (bottom) for 1200 UTC 28 February 1991. 
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Figure 15. Wqf ior 1200 UTC 28 February 1991 and observed precipitation 
(inches, bottom), 0600 UTC 28 February 1991 - 1800 28 February 1991. 
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(iTum, top) and vrO/ (bottom) for 1200 UTC 04 March 1991. 
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Figure 17. Wqf {top) for 1200 UTC 04 March lyyi and observed precipitation 
(inches, bottom), 0600 UTC 04 March 1991 - 1800 UTC 04 March 1991. 
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Figure 18. Unsmoothed (top) and smoothed (bottom) orographically forced 
vertical velocity (cm/s). 0000 UTC 01 March 1991- 
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RE(?TONM 



1. App»'pach 

While the general results are indicative of the overall value of the wqf 
parameter in representing precipitation over a wide region, it is perhaps more important 
to determine the extent to which the correlation can be improved when applied to a 
specific topographic region. In tliis way, wqf may be able to identify specific regions or 
even diainage basins likely to receive heavy precipitation under certain forecast 
conditions. For this reason, the specific perfonnance of wqf in regions of well defined 
topography is examined to determ.ine the ability of wqf to represent precipitation in 
specific regions. In effect, the extent to which synoptically derived low-level winds 
interacting with topography can represent the mesoscale distribution of precipitation 
within mountainous regions will be determined. Little or no improvement of the regional 
coiTelation coefficient over the general case would suggest limited usefulness of the 
synoptic-scale winds in determining orographic forcing. An increase in correlation would, 
however, indicate that a significant component of the precipitation forcing is due to the 
large-scale wind field interacting with regional topography, and that this quantity is more 
easily isolated in regions of well defined topography. Specific correlation coefficients are 
presented as required in each sub.region. 

The regions chosen for this specific purpose all have topography which is 
relatively well represented at 40 km resolution. In each case, regions are defined by the 
particula.r topographic fean.Mes ^ < ' ' Res’ilt<: are 
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presented by showing the obseirved precipitation across a topographic cross-section which 
is chosen based on observational density as well as topographic gradient. Observations 
are then compared to the derived quantities. Topography requiring finer resolution to 
be represented does not fare well, as smoothing tends to destroy the forcing maxima (a 
good example is the Santa Cruz mountains). As much as possible, limitations, as well 
as strengths, corresponding to specific regional correlation will be identified. 

2. Topographic Subregions 

o. Son Bernardino/San Qobriel 

This subregion was particularly hard hit by storms early in the period. 
Figure 20 shows the regional topography to consist of a broad flat area bordered to the 
northeast by the abrupt rise of the San Bernardino range. Figure 21 shows the actual 
observed precipitation accumulations along with the derived quantities, for 0600 UTC 
-1800 UTC 01 March 1991 combined with a representative cross-section of the regional 
topography. This particular cross-section is representative of the general correlation 
within this subregion, and was chosen since it represents the period of maximum rainfall. 

The wqf parameter seemed to do a particularly good job handling the 
Mojave Desert stations throughout the period, while showing only moderate success in 
the upslope regions. This may best be explained by the coarse resolution of the 
topography southwest of the crest. Our technique sees only a broad, continuously sloping 
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Bemardino/San Gabriel subregion (0000 UTC 28 February 1991, top) and Southern 
Sierra Nevada subregion (1200 UTC 01 March 1991, bottom). 
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region, whereas many reportjjrg stations lie in dissected valleys (i.e. Riverside at 301 m 
is assigned an elevation of 550 m by the 40 km topography). Also note (stations 1, 2 and 
3 on Figure 21) the inflated wqf values close to the ocean. This results from the 
smoothing process described earlier. 

With some notable exceptions (i.e., Lytle Creek Ranger Station (LYCR 
#11)), Mv^fdoes a decent job of representing the regional, precipitation showing a 10% 
increase in the correlation coefficient over the state-wide case. Forcing responsible for 
the anomalous observations (i.e., LYCR and CRES) is certainly of too small a scale to 
be accurately accounted for by this technique. Still, the ability of wqf to correcdy 
produce a specific area of heavy precipitation (>1.5" in 12 hours), as well as define the 
lesser amounts observed in the Mojave Desert (stations VICT, BORO, DAGG and MOJA) 
is encouraging. By comparison, NGM 12 hour precipitation totals forecast for this region 
vvere all less than 0.4 inches with no upslope or lee structure in the NGM precipitation 
(not shown). 

Of course, wqf has the benefit of being statistically fitted to actu<d 
observations, where the NGM forecast represents a strictly physical quantity. 
Nonetheless, the statistical fit achieved by u’^is derived objectively, and is intended to 
be applied in all conditions. Still, an application of this method to independent data is 
needed to determine predictive skill. 
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Figure 29- Topography (meters) and station names for San Bemardino/San 
Gabriel subregion. Precipitation amounts (inches) are for 12 hours ending 1800 UTC 
01 March 1991. 
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Figure 2.1. Topography (solid line, meters), precipitation (dashed, inches), 
and observed precipitation (points, inches) for San Bemardino/San Gabriel subregion 
(12 hour period ending 1800 UTC 01 March 1991). 



It should be mentioned that the significant rain events in both this 
subregion and the Santa Barbara subregion which follows ended by 1800 UTC 02 March 
1 991. As a result, the conditions which lead to deterioration of correlation values late in 
the period were not obse.rved in these regions. 

b. Santa Barbara 

The Santa Barbara subregion (Figure 22) is characterized by particularly 
good H'^Z/piecipitatiou '.cnelof'o'' j ’■> -'i, . ,.,.,1 . ... m rouging from 
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-0.52 to -0.85 for the different analysis times. Note that correlation during the post- 
frontal convective period of 1200 UTC 28 February 1991 - 0000 UTC 01 March 1991 
shows less deterioration in this region than in the state-wide case (-0.52 vs. -0.22). 
Similar to the previous subregion, lee stations (NE) are represented highly accurately, 
with the majority of the variance seen in the upslope regions (SW). Again, subgrid-scale 
topography likely produces large variations over small distances. 
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A good example of ajn error caused by insufficient resolution of the 
topography is evident at time 0000 UTC 28 February 1991. Figure 24 shows the wqf 
maximum to be "transposed” to the southeast from the observed maximum, therefore 

litniting any coixeJatio^^ in in this 
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Figure 23. Topography (solid line, meters), precipitation (dashed, inches), 
and observed precipitation (points, inches) for Santa Barbara subregion (12 hour 



period ending 1800 UTC 01 March 1991. 



region (not shown) depicts a narrow northeast trending ridge in the vicinity of Thousand 
Oaks (station LECH), northwest of which is the Ventura Valley (station SEMI). This 
topography is represented as a continuous slope at 40 km resolution, however. As such, 
the strongest orographically forced component of wqf is located over the Ventura Valley 
(near SIMI), whereas in reality this forcing maximum is located between the Ventura 
Valley and Pine Mountain (station PINM). As the southeasterly winds become more 

southerlv after this tw.tiH 
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Figure 24. Wqf and observing stations with 12 hour precipitation totals ending 
0600 UTC 28 February 1991 (inches) for the Santa Barbara subregion . 

favorably, as the large-scale winds interact with topography that is better represented at 

the 40 km resolution. 

Overestimation of precipitation at the coast (Figure 23, stations 1-5) does 

not seem to be as p*'e'''a1e’’' thi- i'.. . .-.i- ;> oi 

slope at the )and-ocea't Si'ectfical^v. '''■''''.'•Hint; P” topographic u‘- 
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field still produces upward orographic forcing well into oceanic regions, the rapid increase 
in elevation onshore produces orographically enhanced observed precipitation amounts 
even near the ocean. Moreover, the apparent high coastal correlation is more likely a 
coincidence as opposed to an increase in skill over, for example, the San Bernardino 
region. As with the previous sub-region, no significant precipitation was observed after 
1800 UTC 02 March 1991. 

C- Southern Sierra Nerada 

The most notable feature seen in this region (Figure 25), commonly 
known as the High Sierra, is the lack of observations. Still we can make at least one 
important observation: 

- Wqf generally overestimates precipitation for these inland locations. 

This trend is evident throughout the study period and is represented in 
Figure 26. We can generally note that wqf responsible for a specific rain amount along 
the coast (i.e., Santa Barbara) represents a somewhat lesser amount at more inland 
locations. This feature is almost certainly a function of the available moisture, such that 
inland locations require a greater wqf to maintain saturation than is required near the 
coast. No effort is made to quantify this relationship, however it must be considered 
when using wqf based on inter-regional statistics (similar to adjusting Model Output 
Statistic (MOS) precipitation forecasts for a specific location). 
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Figure 25. Topography (meters) and station names for Southern Sierra Nevada 
subregion with 12 hour precipitation totals (inches) ending 1800 UTC 01 March 



1991. 



Note that Figure 26 shows u’^/and observed precipitation both before and 

ahc l200t_rrrp2.Ma>''l' l '''>l v Wh ’ -'i : <>( 0,72 nn<i - 

0.70 respecttvclv. Thonj^b in tli'- Putf'i . .1 'I'- t;'--'- <•* I <>t '• tH’ll- 
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Figure 26. Same as in Figure 23 except for Southern Sierra Nevada subregion 
(12 hour period ending 1800 UTC 01 March 1991, top, and 0600 UTC 04 March 
1991, bottom). 
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orographic effects are more evklent, especially considering the rough equivalency of 
precipitation totals at the San Joaquin Valley station Corcoran (CORC) and amounts 
recorded higher in the Sierra. TTre lovv slope value (-0.004) of the linear regression 
supports this assignment. 

d. Northern Sierra Nevada 

This subregion (Figure 27) depicts the gradually sloping region of the 
Siena Nevada between Sacramento and the Lake Tahoe basin. The chosen cross-section 
(Figure 28) coincides roughly with the American River drainage, which is represented 
particula.rly well by the 40 km topographic resolution. As a result, we see very good 
M'<5';7precipitation conelation in this sub-region prior to 1200 UTC 02 March 1991, with 
correlation coefficients between -0.70 and -0.80. As before, we note the overestimation 
of prec ip.it ation in the Sacramento valley (WHEA), as well as a general tendency to 
underestimate precipitation (slope values average -0.006 in this subregion). The generally 
good fit seen in the finst period at the summit and lee stations of Soda Springs (SODA) 
and Truckee (TRUC) respectively represent the ability of this technique to account for lee 
"spillover", due to the statistical correlation producing a positive y-intercept. 

Perhaps the best example of subgrid-scale processes eiffecting the 
distribution of precipitation is seen at time 1200 UTC 04 March 1991 (Figure 27). Note 
the extraordinarily high precipitation observations over Dormer Summit, including 
Tmckee, as compared to precipitation estimated by ^vqf (regional correlation coefficient 
OIkI Slop^ values PJe 'vvl -O (V' l ,*lv 11*- • • n.naiv,^i ; not 

resolving the win4 m<^M5tnr^ f ; r 1 1 i thijtif’vii. This 
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extimple brings into question u'c//’ s level of usefiibiess given this type of synoptic 
situation. 

Of particular note concerning this region is the ability to compare results 
obtained here to those obtained over a similar region by Colton’s (1976) 4.3 km model. 
Though Colton’s (1976) use of an hourly precipitation rate prohibits direct comparison 
to the 12 hour rate presented here, we can make some qualitative assessments. Other than 
Colton (1976) showing considerably greater detail within the upslope region, including 
the secondary maximum produced by Mt. Rose rising to the east of Lake Tahoe, wqf 
produces a similar enhancement pattern, prior to 1200 UTC 02 March 1991, in terms of 
the location of maxunum as well as tire intensity. Of further note is that even at the 4.3 
km grid scale, Colton’s (1976) calculations at lee stations generally underestimate 
precipitation, even with a "spiUcver" temi iticluded. However, Colton’s (1976) fine mesh 
resolution allows a better representation of convective forcing. In these types of 
situations, wqf is far inferior to Colton’s (1976) model. 
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Figure 27. Topography (meters) and station names for Northern Sierra Nevada 
subregion with 12 hour precipitation totals (inches) ending 1800 UTC 04 March 
1991. 



e. Shast<x/Siskiyoff 

The final sub-region (Figure 29) contajjis not oaly the most complicated 
topogrophv, but U' tli^ witi' ), m-Nm; to pi pripitntion. 

Figure ?0 s^ows rcp.Tese^uc.t.^ve ^’v/Z/piecipiieUi/'i’^ i.-t' • I zno L^C 
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Figure 2S- Same as in Figure 23 except for Northern Sierra Nevada subregion 
(12 hour accumulations ending 0600 UTC 02 March 1991 , top, and 1800 UTC 04 
March 1991, bottom). 
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02 March 1991 with correspondmg correlation coefficients of -0.68 and -0.46 respectively. 
Also, slope values for this subregion average -0.012, which, being relatively high, is 
consistent with its more coastal location. Despite the obvious smoothing of topographic 
features at this resolution, the overall wqf representation of precipitation is quite 
encouraging, in particular the secondary maximum seen between Happy Camp (HAPP) 
and Elk Valley (ELKV). Even more interesting is the less obvious loss of correlation late 
in the period. This may be due to tlris subregion’s northernmost location. Here, 
synoptically derived subsidence is less than that at more southern locations at this time, 
allowing wqf to "create" upward motion by offsetting the synoptic -scale subsidence. 
Clearly, the improved correlation is only coincidental since convective processes are not 
accounted for in any region. Another possibility is more consistent forcing. The 
proximity of this subregion to the surface low pressure systems late in the period may 
have resulted in a more consistent low-level southerly flow. In this situation, though 
convection is still unaccounted for, the synoptic- scale winds are probably acceptable 
representations of the low-level flow in this subregion. This, combined with the strong 
tendency for the north end of the Sacramento Valley to channel winds with a southerly 
component may have allowed wqfio reach a level at least of the same order of magnitude 
as the convective forcing. As such, this period’s correlation is likely more an effect of 
the particular local situation interacting with the synoptic scale, rather than wqf being 
more accurate in this subregion 
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Figure 29. Topography (meters) and station names for Shasta-Siskyou subregion, 
with 12 hour precipitation totals (inches) ending 1800 UTC 01 March 1991. 
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Figure 30- Same as in Figure 23 except for Shasta-Siskyou subregion (12 hour 
accumulations ending 1-SOO UTC 01 March 1991, top. and 1800 UTC 04 March 
1991. bottom). 
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f. Non -mountainous Regions 



Though this study is aimed at representing mesoscale precipitation within 
mountainous regions, it is also hoped that including the effects of local topography will 
improve distribution in low lying regions, or regions in which topography is sub-grid 
scale such as the Santa Cruz Mountains (Figure 31). Without going into great detail, we 
can qualitatively assess the success of wqf as weaker where topography is not as well 
defined. As before, regions of less dramatic topography are most susceptible to sub- 
synoptic processes dominating the local distribution of precipitation. This is especially 
evident vvhen surface forcing is not represented well by the synoptic analysis (i.e. upper 
level features and convection). 

As far as mountainous regions which are poorly represented at 40 km 
resolution, there is sufficient evidence to suggest improving topographic resolution in 
order to identify wqf forcing in these regions. Again we refer to the precipitation 
distribution within the Santa Cruz mountains (Figure 32). Present resolution does not see 
the steep, narrow nature of this range, but observations support strong local orographic 
influence. Increasing topographic resolution may provide sufficiently large quantities of 
wqf to "survive" the smoothing process realizing that this representation wUl only be of 
use if the synoptic analysis is at least close to reality, which is certainly far from 
guaranteed. 
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Figure 31. Topography (meters) and station names in the vicinity of the Santa 
Cruz Mountains with 12 hour precipitation totals (inches) ending 1800 UTC 28 
February 1991. 

3- Comparison to NGM Precipitation Forecasts 

Earlier we established the rough equivalency of the NGM precipitation forecast 
with the term used here to represent synoptic-scale forcing, wOf. We have seen examples 
of the improvement the topographic forcing term has on the wOf, field in representing 
observed precipitation during the first period. This is essentially the same level of 

improvenient vill hf-ei' ."Iv'v. ''\ >m ih*- i i' -r ’ ' .< v. i»ii !i- ■ • i-iK'n allowing 

for the NGM convectiv*;' parani^ipi (this pa,’a.o''''0;' • jn-t n : . n i'\- hnilt 



66 




1800 28 February 1991 (inches) in vicinity of Santa Cruz Mountains. 



into wqf). Note that while the NGM and wO/are able to define general precipitation areas 
reasonably well, there is virtually no structure other than that associated with the synoptic- 
scale forcing maxima (usually coincident with 700 mb vertical velocity), and therefore 
weak correlation with observations. 

Late in the period. NGM represents precipitation with only slightly less overall 
accuracy than wqf. Recall that for these times, the NGM vertical velocity forecast 
predicted upward motion in areas which verified to be subsident. Still, its convective 
parameter may be • • i--. - 
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VI. CONCLUSIONS 



A. APPLICATION TO FORECASTING 

J.t is -well known that precipitation amounts can var>' dramatically over a small area. 
Inclusion of topography into fine-mesh models has shown favorable representation of this 
variability (Colton 1976, Collier 1975). Unfortunately, the data network needed to run 
these models is not generally available to the forecaster, who can only combine local 
knowledge witli the synoptic-scale products he receives. In this thesis we have attempted 
to provide a quantitative esti.nrate of mesoscale precipitation distributions over 
mountainous terrain, without requiring mesoscale infonnation, other than topography. 
Conclusions can be summarized by the following points: 

* Given that synoptically mduced forcing is well represented in the analyzed or 
forecast fields, wqf can predict a substantial degree of the mesoscale precipitation 
distribution. 

* Most of the discrepancies between wqf and the observations occur within upslope 
regions (as opposed to lee stations) due to unresolved small-scale topographic 
forcing. 

- Correlation is greatly improved when considering well defined topographic regions. 

* Wqf shows significant improvement over the NGM precipitation forecast fields. 

Coincident with our conclusion that wqf is a good indicator of precipitation under 
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Therefore, to be of use to the forecaster, he/she must fust determine the applicability of 
MYjrf based on the synoptic situation. There is sufficient evidence to suggest w^may be 
of some value even when no synoptic -sc ale forcing is indicated. In this case, the 
forecaster must weight the wqf field accordingly, but be aware that areas of consistently 
strong flow may still be represented by wq/' precipitation maxima. As such, assuming wqf 
is produced as a forecast field, the forecaster must identify its credibility based on the 
following parameters: 

* Presence of a well-defined thickness ridge ahead of a front. 

- The potential for unstable precipitation. 

If these conditions can be correctly evaluated, wqf will be sufficiently useful to the 
forecaster in predicting locally heavy precipitation to warrant its calculation by regional 
forecast offices. 

B. RECOMMENDATIONS 

The most immediate recommendation would be to incorporate actual NGM forecast 
fields into the wqf calculation. This would provide a basis for verifying the wqf 
precipitation against the actual NGM forecast precipitation amounts. Such an application 
will be able to more definitively evaluate wqf as a forecast parameter. 

Secondly, there may be some benefit to computing wqfovtr. even finer topographic 
resolution (i.e.. 20 c in i ,>,> -jiiv * <'i. (•,. i... -- n;. i< — ••• i-'l-iii-Mi will be 

prevalent regardless 0^^ *h.e i.ono^' anhic te."olu'i''” '•’''■ii a” ii'C!-'> •; 'es'^l'.'iinn iriov 
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improve wqf representation, of precipitation by pinpointing areas of maximum 
ejtliancement and by reflecting the orographic influence of steep topography with less 
horizontal extent. There will, however, most certainly be a point at which no increase in 
wqf accuracy is seen with increasing topographic resolution. Furthermore, objective 
restrictions potentially can be applied to the calculation of wqf in flat regions which are 
adjacent to significant topographic barriers. The intent here is to limit the forcing to the 
actual topographic slope rather than distributing it over a larger region as done by the 
crment smoothing process. 

As a final note, I wish to emphasize that no new ground has been broken here as 
far as our understanding of the precipitation process. I have simply adapted seeder-feeder 
theory of orographic ally enliajtced precipitation to a format which only requires synoptic 
scale information. Though the results are far from ideal, I think the technique exhibited 
sufficient success to warrant its consideration by the forecaster. 
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